NASA Technical Memorandum 83794 


Ultrasonic Velocity Measurement 
Using Phase-Slope and Cross- 
Correlation Methods 


David R. Hull, Harold E. Kautz, 
and Alex Vary 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

1984 Spring Conference of the American 
Society for Nondestructive Testing 
Denver, Colorado, May 21-24, 1984 





E- 2290 


ULTRASONIC VELOCITY MEASUREMENT USING PHASE SLOPE AND 


CROSS CORRELATION ME I HODS 


David R. Hull, Harold E. Kautz, and Alex Vary 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 441 3b 


SUMMARY 

Computer Implemented phase slope and cross- correlation methods are In- 
troduced for measuring time delays between pairs of broadband ultrasonic 
pulse echo signals for determining velocity In engineering materials. The 
phase slope and cross-correlation methods are compared with the overlap method 
which Is currently In wide use. Comparison of digital versions of the three 
methods shows similar results for most materials having low ultrasonic attenu- 
ation. However, the cross-correlation method Is preferred for highly attenua- 
ting materials. An analytical basis for the cross-correlation method Is 
presented. Examples are given for the three methods Investigated to measure 
velocity in representative materials in the megahertz range. 


1NTR0DUC I ION 

Ultrasonic velocity measurements are widely used to determine properties 
and states of materials. In the case of engineering solids measurements of 
ultrasonic wave propagation velocities are routinely used to determine elastic 
constants (refs. 1 to 5). There has been an Increasing use of ultrasonic 
velocity measurements for nondestructive characterization of material micro- 
structures and mechanical properties (refs. 6 to 10). Therefore, It Is Impor- 
tant to have appropriate practical methods for making velocity measurements on 
a variety of material samples. 

A number of approaches to the measurement of ultrasonic velocity In solids 
have been described (refs. 11 and 1?). Included among the approaches there are 
those that use single frequency continuous waves, tone bursts and broadband 
pulse echo waveforms. In the latter case, the basic problem Is to determine 
the time delay between two successive echoes, each of which consists of only a 
few oscillations. The pulse-echo approach’ using a single transducer In contact 
with the material sample has many practical uses (refs. 13 and 14). Therefore, 
It will be treated exclusively In this report. 

In many Instances It Is fairly easy to find some feature in the pulse-echo 
waveforms (e.g., the principal oscillation) that can be used as a basis for 
measuring the time between echoes. Ihe "pulse-echo overlap" method (ref. 11) 

Is commonly used to make this type of measurement. The method can be Imple- 
mented by analog or digital Instrumentation (ref. 12). In either case, the 
echoes must have similar waveforms so that corresponding features can be read- 
ily Identified and brought into coincidence (l.e., overlapped) either by analog 
or digital manipulations. Since this condition Is not always met, alternative 
methods arc needed. 



The purpose of this report Is to deal with those cases where pulse-echo 
signals are weak or distorted by attenuation and other factors that render them 
unsuitable for the overlap approach. This report describes ultrasonic velocity 
measurement methods based on computer digitization of broadband pulse-echo 
waveforms. Three digital methods are compared: (1) overlap, (2) phase-slope, 

and (3) cross-correlation. The analytical basis for the cross-correlation 
method Is emphasized. Specific examples are given for observed signals In a 
frequency range from 5 to 50 MHz for metal and composite samples. 


EQUIPMENT 

Transducer 

A typical transducer configuration for normal Incidence, contact, pulse- 
echo ultrasonics appears In figure 1(a). In the transducer a piezoelectric 
crystal Is excited by a high voltage drive pulse to produce a broadband longi- 
tudinal ultrasonic wave. This wave propagates through the quartz buffer to 
the front surface where a portion of the wave's energy Is transmitted Into the 
specimen and a portion Is reflected back to the piezoelectric crystal. The 
portion of the ultrasonic wave reflected from the front surface Is the first 
echo, FS (fig. 1(b)). The wave transmitted Into the specimen Is reflected by 
the back surface to the front surface where a portion Is reflected back Into 
the specimen and a portion Is transmitted through the quartz buffer to the 
piezoelectric crystal. This wave Is shown as echo B1 In figure 1(b). Finally 
a third ultrasonic wave, echo B2, occurring from a second back surface reflec- 
tion Is received at the plezoelctrlc crystal as shown In figure 1(b). 


Instrumentation 

Block diagrams of Instrumentation used for this report appears In 
figures 2(a) and (b). The major components are: (1) ultrasonic pulser- 

recelver, (2) digitizing oscilloscope, and (3) minicomputer. Auxiliary com- 
ponents Include a digital CRT monitor, real-time video monitor, digital time 
delay, and hard copier. The minicomputer operates under TEK SPS BASIC software 
(ref. 15) and various system components are Interfaced with the IEEE 488 6PIB 
(General Purpose Interface Bus). The ultrasonic pulser-recelver Is an analog 
device used to excite the transducer and to receive resultant echoes. There- 
after, signal acquisition, processing, and analysis are done digitally. Meth- 
ods for digital processing described herein depend on use of the discrete 
Fourier transform algorithm. It Is Implemented by an FFT (fast Fourier trans- 
form) software package which Is part of the system of figure 2. 


PROCEDURE 

General 

In figure 1(b) time, T, between echoes B1 and B2 Is that required for 
ultrasonic pulses to travel twice the thickness, s, of the sample. The time 
between echoes FS and B1 Is also T. Here echo FS will be ignored because It 
Is preferable to use echoes from within the specimen. The procedure for calcu- 
lating velocity, v = 2s/T, Is based on analysis of the two back surface echoes 
B1 and B2. All examples given are based on longitudinal wave propagation. 
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The Initial step In the procedure for signal analysis requires digitiza- 
tion of echoes B1 and B2. Two typical cases can arise: In the first case, 

both echoes are displayed simultaneously on one CRT trace (as In fig. 2(a)). 

In the second case, the echoes are Individually “windowed" In two separate 
traces (as In fig. 2(b)). These "windows" are time gates used to Isolate wave- 
forms of Interest. In digitized form they usually consist of 512-element 
arrays with 9 bit vertical resolution. The Individually "windowed" case Is 
used when there Is a long time Interval between echoes. Then, any attempt to 
compress two echoes Into a single CRT trace can significantly reduce the amount 
of usable detail In the Individual echoes. The signals could also be rendered 
useless due to aliasing, according to the Nyqulst criterion for digitization 
(ref. 16). 

When both echoes B1 and B2 cannot be simultaneously windowed, a digital 
delay device Is used. It Introduces a known time delay between the beginnings 
of the two separate windows containing the waveforms for echoes B1 and 82. The 
time delay device Is calibrated to an accuracy of ±1 ns and has a precision of 
±0.1 ns. As Indicated In figure 2(b), this delay Is triggered by synchroniza- 
tion signals from the pul ser-recel ver . After B1 Is windowed and digitized, a 
predetermined delay Is generated to put echo B2 Into the digitizing window. 

The process Is shown In figure 3 for a case where the time, T, between 
echoes B1 and B2 Is calculated from T = W * (T2 - Tl), where, W Is the time 
delay between windows and times Tl and T2 are In-window delays measured 
with respect to corresponding major oscillations of echoes B1 and B2, respec- 
tively. If both echoes can be Included In the same window, W = 0 and the 
expression reduces to T = (T2 - Tl). 


Overlap Method 

Figure 4 shows the result of overlapping the echoes B1 and B2 from 
figure 3. The process Involves a digital search for the peak voltage and the 
array element that contains It for each echo's major oscillation. Since, ex- 
cept for amplitude, the echoes have similar waveforms the peak value criterion 
Is adequate to overlap the echoes. 

The overlap method depends on having a pair of echoes, B1 and B2, that 
exhibit similar waveforms with corresponding features. But, even for the case 
Illustrated In figure 4, It Is apparent that the overlap Is Imperfect, e.g., 
the zero crossings are not all coincident. Effects of wave distortions due to 
noise, dispersion, and other factors that operate on successive echoes are 
present. The ability to deal with these distortion effects Is facilitated by 
operating In the frequency rather than the time domain. It will be seen that 
frequency domain methods (phase-slope and cross- correlation) have advantages 
over the time domain method. 


Phase- Slope Method 

With the phase slope method time between echoes Is found by use of phase 
spectra of echo waveforms. After the echoes are digitized, a Fourier transform 
of each Is obtained by a discrete F FT algorithm. The amplitude and continuous 
phase spectra for a pair of typical echoes are Illustrated In figure 5. 
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After Fourier transformation, both the amplitude and phase spectra are 
used to define a central zone within the frequency domain. For example, this 
zone may consist of only a narrow range near the center frequency or a fre- 
quency range for which the amplitude exceeds some fraction of the peak value. 
And/or, the zone may consist only of the frequency range for which the phase 
spectrum Is linear. These restrictions eliminate the low and high frequency 
extremes where the slgnal-to-nolse ratio Is low. 

The group velocity Is given by U(f) = 2w(2s)/((de/df ) - W) where f 
Is frequency, 0 Is phase angle In radians and W Is the time delay between 
echo windows (ref. 17). In window time delays for each echo can be calculated 
/ by setting W = 0 and solving for T(f) = (de/df )/( 2ir) . In the present 

case the phase spectra are linear functions of frequency In the central zone, 
their slopes are constants, M = de/df, and the In-window time delays are 
T1 = Ml /( 2ir) for echo B1 and T2 = M2/(2«) for echo B2. The total time delay 
T . W * (T2 - Tl). 


Cross-Correlation Method 

The digital cross-correlation method eliminates the need for the somewhat 
arbitrary criteria, e.g., peak value and zone for phase slope, applied In the 
two previously described methods. Unlike the overlap or phase-slope methods, 
cross-correlation does not require explicit criteria for accepting or rejecting 
specific features In echoes affected by distortion or low slgnal-to-nolse 
ratios. 

The method Is Illustrated In figure 6 which shows the normalized cross- 
correlation function for three cases of echo positions. The cross-correlation 
function of two time domain signals is obtained by conjugate multiplication 
using their Fourier transforms In the frequency domain and retransformation 
back to the time domain (ref. 16). The cross-correlation function possesses a 
maximum In the time domain. The displacement of this maximum relative to a 
zero reference gives the time Interval C, which for the Ideal case should 
equal T2 - Tl as measured by the digital overlap method. Cross-correlation 
gives this quantity whether the echoes appear In the same or separate windows. 
If the echoes B1 and B2 are separated windowed, then T = W + C. 

For the relatively simple and undistorted echoes of figure 6 cross- 
correlation gives results similar to those obtained by the two previous meth- 
ods. The validity of the cross-correlation method depends on the fact that 
the displacement of the maximum of the cross-correlation function equals the 
time delay between the two successive echoes. This Is supported by the previ- 
ous examples and the theory and results presented hereinafter. 


THEORY 

The theoretical analysis proceeds according to the following logic. 

First, It Is noted that the center value of an auto-correlation function Is 
always a maximum. Then, It Is shown that the cross-correlation function can be 
rewritten as an auto-correlation function, this rewritten form has a maximum 
displaced from Its center and that the displacement corresponds to the time 
delay sought 
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The cross-correlation function R X y( T ) defined In reference 18 as, 


V T) 


11m 

T-x*> 


( 1 / 1 ) 


/ 


X(t)Y(t + x)dt 


( 1 ) 


where, X(t) and Y(t) are arbitrary waveforms, where t Is time and t Is 
a delay Interval. If Y(t) Is replaced by X(t), we have the auto-correlation 
function, 


R vy (t) = 11m (1/T) 
X * T -»00 



X(t)X(t * -t)dt 


( 2 ) 


It Is assumed that the data of Interest are stationary, l.e., that the 
time averages of equations (1) and (2) exist. It Is shown In reference 18 
that R xx (0) - R(t) xx f° r T - That Is, the center value of the auto- 
correlation function at t = 0 Is always a maximum. 


In an analysis of propagation paths (refs. 16 and 18) (as those associated 
with echoes B1 and B2) a function can be defined. This function operates 
along the kth propagation path on the Input signal X(t). The resultant output 
signal Y(t) Is, 


y(t) = lH k X(t)fA(t) (3) 

where A(t) represents accumulated random noise. Equation (3) applies to 
pulse echoes (as In fig. 1(b)) If the Input function X(t) Is taken as the 
first echo B1 . Subsequent echoes such as B2 are functions of the propagation 
function Hfc. By excluding higher order echoes, equation (3) reduces to, 

Y(t) = HX( t) f A(t) (4) 

In the case Illustrated In figure 1(b) the function H operates on the 
Input signal X(t), or echo B1 , In the following ways before It yields Y(t), 
or echo B2. First, It attenuates echo B2 relative to echo B1 . Second, It may 
Invert echo B2 relative to echo B1 (e.g., reflection at a free back surface, 
l.e., multiplication by -1) followed by a second Inversion at the buffer-sample 
Interface, reinverting B2 to give results as In figure 1(b)). Third, It delays 
echo B2 relative to echo B1 by the round trip delay time In the sample. The 
delay time, T, Is related to the propagation velocity, v, and sample thickness, 
s, by T = 2s/v. The effect of H operating on X(t) can be expressed as: 

HX( t) = ±DB(f )X(t) (5) 

The ± sign on the right side of equation (5) corresponds to Inversions 
of echo 82 relative to echo B1 and depends on reflection coefficients at the 
buffer/sample Interface. The time delay 0 Is an operator changing t to 
t - 2s/v. B Is the attenuation operator which Is, In general a function of 
frequency f. B must operate In the frequency domain as shown: 
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BX( t) = F" 1 [BF(X( t) ) ] 


(5a) 


Here F Is the Fourier transform and F -1 Is the Inverse Fourier transform, 
defined so that: 


F" 1 [F(X(t)>] (5b) 

If, however, B does not distort the shape of X(t) too much, then we can 
approximate It by a constant b and write; 

Y(t) = HX( t) = ±bDX( t) 

= ±bX(t - 2s/v) (6) 

According to reference 15, the discrete cross-correlation function for a digi- 
tized finite waveform can be written as, 

N-n-1 

R vv (n) = (l/N) Z X(t)Y(n + t) (7) 

xy t=0 

where N equals the number of data sampling Intervals corresponding to the 
delay time T. The quantity n ranges from -(N - 1) to +(N - 1). Combining 
equations (6) and (7) for the cross- correlation of echoes B1 and B2 gives, 

R (n) = (±)b/N£x(t)X[n + t - (2s/v)] (8) 

* J 

In reference 16 It Is shown that equation (8) Is equivalent to the auto- 
correlation of X(t) with the correlation array points shifted to the right a 
distance 2s/v. Thus the central maximum, or minimum, Is shifted to the point 
a distance 2s/v right of center, In the case of a time lag (delay) from the 
first to second echo (as In fig. 6(b)). 


EXAMPLE 

The previous examples of digital methodology (figs. 3 to 6) were based on 
a metal sample (250-grade maraglng steel). The velocities, measured with a 50 
MHz transducer, are compared In table I for each of the three digital methods 
used. As expected, the results agree quite closely for the relatively clean, 
uncomplicated waveforms Involved. 

A further experimental example Involving noisy, complicated waveforms Is 
Illustrated In figure 7. In this case both echoes, B.l and B2, occur In the 
same window. The signals were obtained from a graphite/epoxy composite lami- 
nate using a 5 MHz transducer. Due to the nature of the material, there Is 
considerable attenuation, distortion, and noise. In addition, echo B2 Is In- 
verted relative to echo B1 because of the sign of the reflection coefficient 
at the buffer/sample Interface. 

The example In figure 7 demonstrates the facility of the cross-correlation 
method. Figure 7(a) shows echoes B1 and B2 and figure 7(b) shows the corre- 
sponding cross-correlation function. In figure 7(b) the maximum Is displaced 
from 0 by 0.750 ± 0.002 ys. It Is negative due to the Inversion of echo B2. 
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figure 7(c) exhibits highly distorted and noisy echoes obtained by 
windowing the third and fourth back surface echoes, B3 and B4, without removing 
the specimen from the transducer. The. time between this pair should be Identi- 
cal to that between echoes B1 and B2.. Therefore, the result of applying the 
cross-correlation method In figure 7(d) can be compared directly with the pre- 
viously determined value for time T. For figure 7(d) T = 0.752 ± 0.002 ps 
with the difference between the two within the timing error corresponding to 
the data sampling Interval. However, the digital overlap method gave 
T = 0.748 ± 0.002 ps for echoes B1 and B2 and T = 0.766 ± 0.002 ps for echoes 
B3 and B4. The computer program used to the results In figure 7 Is listed In 
the appendix. 

Even In the case of clean, uncomplicated echoes, all three methods are 
subject to a "data sampling" error peculiar to the digitization process. This 
Is due to the finite number of discrete data points that the digitization proc- 
ess assigns to the waveforms collected. A random error tends to arise from a 
variety of sources, e.g., small jittering effects In the oscilloscope sweeps. 

To estimate this Inherent, random error, a series of 30 waveform pairs were 
digitized. In sequence, at the same location on a material sample. The results 
appear In table I which compares the three digital methods for a metal sample 
and a composite sample. Typical errors In the time delay are given In table I 
for the particular material samples, frequencies, and signal processing methods 
used In this study. Note that the time delay errors are compared for a fixed 
thickness of each material. It Is seen In table 1 that the cross-correlation 
method gives the lowest random time delay errors for both samples. The phase- 
slope method on the composite sample produced the highest random time delay 
error due to the low slgnal-to-nolse ratio of the echoes. 


DISCUSSION 

The digital overlap method Is an effective way to automate the analog 
overlap method. The overlap method requires Input of Information concerning 
echo Inversions, peak selection, etc. Once overlap criteria are established, 
the digital time domain method simulates a manual overlap based on visual In- 
spection of analog waveforms. However, even digital signal manipulation of 
time domain signals can be problematic If the waveforms are subject to noise, 
attenuation, and similar factors. 

The frequency domain phase- slope method eliminates problems encountered 
In the time domain, e.g., the need to account for echo Inversions. In addi- 
tion, It provides convenient criteria for selecting an appropriate frequency 
range In cases where the major portion of the phase spectra of the echoes are 
mutually linear. Generally, for nonlinear dispersive cases, the phase-slope 
method can determine group velocities as functions of frequency (ref. 17). 
However, If slgnal-to-nolse ratio Is low, as was the case for the composite 
sample In table I, poor results are obtained. 

The advantages of the cross-correlation method are apparent when the 
slgnal-to-nolse ratio Is low and/or random noise Is superimposed on the echoes. 
One of the demonstrated properties of the cross- correlation function Is that It 
Is (statistically) weighted by dominant frequencies common to the waveforms 
being correlated (refs. 16 and 18). Therefore, the cross-correlation function 
returns a group velocity within the frequency bandwidth of the signals ana- 
lyzed. With the cross-correlation method, effects of random noise, extraneous 
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signals, etc. are minimized:' The example of figures 7(c) and (d) Illustrates 
the utility of the cross-correlation method for highly attenuating and noisy 
materials . 


CONCLUSION 

Three computer Implemented digital methods for measuring ultrasonic veloc- 
ity In solids are described: the overlap, phase-slope, and cross-correlation 

methods. It Is evident from the examples given that the cross-correlation 
method Is the most versatile of the three digital methods for measuring time 
delays between broadband ultrasonic pulse-echo signals. All three methods 
produce comparable results when applied to pairs of clean, undistorted echoes. 
The cross-correlation method yields the most reliable and reproducible results 
In the case of highly attenuating materials that return very noisy, distorted 
echoes. A discussion of theory and relevant properties of the cross- 
correlation function Indicates the basis for the versatility and utility of 
the cross-correlation method. 
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APPENDIX 


PROGRAM LISTING 


The program listing given here Is written In TEK SPS BASIC for Implementa- 
tion with a Tektronix 7854 digital oscilloscope and 4052 minicomputer (ref. 7). 
The Instrument configuration Is that shown In figure 2(a) using a Panametrlcs 
5052PR pulser-recel ver and 5 MHz transducer. This Information Is given only 
for Illustrative purposes and does not represent an endorsement of the commer- 
cial equipment or software named. 


100 REM CALCULATION OF UELOCITY BY USE OF CROSS-CORRELATION 
103 REM 

110 REM INSTRUMENTATION: TEKTRONIX 7834 DIGITAL OSCILLOSCOPE 

113 REM 4032 COMPUTER 

120 REM ROM PACKS R07 SIGNAL PROCESSING No. 1 

125 REM R88 SIGNAL PROCESSING No. 2 <FFT> 

130 REM PAHAMETRICS 5052PR PULSER RECEIUER 

135 REM QUART2 BUFFERED TRANSDUCER 

140 REM 

145 REM 

150 REM I INITIAL INPUTS I 

155 REM — — 

160 PRIHT -LEHTER DATE EDO- MHM-YY3 : "| 

165 INPUT D* 

170 PRINT "ENTER TRANSDUCER FREQUENCY : "» 

175 INPUT T* 

180 PRINT "ENTER SPECIMEN ID. : "1 
185 INPUT S* 

190 PRINT ‘ENTER THICKHESS IN cn - "f 
195 IHPUT T 

200 REH 

205 REM I COLLECTION OF WAVEFORM FROM 7854 I 

210 REM 

219 OH SRQ THEN 905 

220 PRIHT "LPOSITIOH ULTRASONIC WAVEFORM ON 7854 WITH" 

225 PRINT "B1 left OF CENTER AND B2 right OF CENTER." 

230 PRINT "PRESS <RETURN> TO COLLECT WAVEFORM f!" 

235 INPUT Z% 

240 V«-"EXEC OFF I 0 2 4 >P/H STORED 1 0 AUG 0 WFM MEAN - 0 WFM SENDX’ 

245 PRINT 01:Vt 

250 IHPUT 0i:AI,A2, A3,A4,A5 

233 DELETE A 

260 DIM ACAD 

263 INPUT 91: A 

270 A«A*A3 

273 A-A+A4 

280 U$-"U" 

285 Ht«*uS" 

290 PRINT 81 .’"SCOPE" 

295 REH 

308 REM I GO PLOT ULTRASONIC WAUEFORM I 

305 REM 

310 PAGE 

315 VIEWPORT 60,120,60,95 
320 GOSUB 930 

325 REM - 

330 REM I BREAK ARRAY IN HALF TO PREPARE FOR CROSS-CORRELATION I 

335 REM 

340 Al*Al/2 

345 DELETE X,Y 

350 DIM X<A1),Y(AI) 

355 FOR J»1 TO A1 
360 X< J) a A< J) 

365 Y< J>=A< J*A1 ) 

370 NEXT J 

375 REM — — 

380 REM I FIND MAX AND MIN OF TIME DOHAIN I 

385 REM 

390 CALL "MAX", X, Ml, Cl 
395 CALL "MAX", Y,M2,C2 
400 CALL "MIN" , Y, M3, C'3 
405 C2 3 C2*Al 
410 13-C3+AI 
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4 1 3 Pun 

420 REM I PERFORM FOURIER TRANSFORM OH EACH HALF FOR PHASE SLOPE I 

425 REM 

430 DELETE 01 , 02, PI ,P2,M 

435 DIM Bl<Al>,B2<Al>,PHAi/2M>,P2<At/2+!>,H<Al'2*l> 

440 B1*X 
445 B2-Y 

459 CALL "FFT",8l 
455 CALL "FFT",B2 

460 CALL "POLAR", 81, M, Pi, 0 
465 CALL "POLAR", B2,M,P2, 0 
470 CALL "DIF3",Pi,Pl 

475 CALL "DIF3",P2,P2 
480 CALL "MAX" , M, P4, P3 
485 P4-1/-<A3*A1) 

490 IF PKP3X0 THEN 500 
495 P1-P1-2FPI 
500 IF P2<P3><0 THEN 525 
505 P2-P2-2FPI 

510 REM - - - 

515 REM I PERFORM CROSS-CORRELATIOM, FIND MAX 1 HIM I 

528 REM - 

525 CALL "CQRR",X, Y,A 
530 DELETE X,Y 
535 CALL “MAX" , A, M8, C8 
540 CALL "MIN" , A, M9,C9 
545 IF -N9<M8 THEM 555 
550 M8--M9 
555 A-A/M8 
560 MF*"" 

565 CALL "MAX", A,M4,C4 
570 CALL "MIH", A,M5,C5 
575 IF-'" 

580 REM . 

585 REM I MAKE DECISION OH IMMERSION fc SET APPROPRIATE UARIABLES I 

590 REM 

595 IF M4>ABS<M5> THEN 620 
600 IF-"<IHMER9I0M>" 

60S M4-M5 

610 C4-C5 

615 C2-C3 

620 HOME Cl, MI 

625 DRAM C1,M1*1.15 

630 MOME C2, M2 

635 DRAM C2.M2FI. 15 

640 REM 

645 REM I GO PLOT CROSS-CORRELATION HAMEFORM i DRAM DELAY MARK I 

650 REM 

655 AI-A1F2 

660 UIEMPORT 60,120,10,40 

665 GOSUB 930 

670 Al-AI/2 

675 MOME C4.M4F1.1 

680 DRAM C4.M4F1.5 

685 PRINT USING "1A, 1D.4D,3A": "<",C4*A3*1000000, " uS" 

690 Tl~(C?-Cl>*A3tl000000 

695 T2-<AUA3*<Pi<P3)-P2<P3)>/P4-'<2tPI)>* 1000000 

708 T3-C4FA3F 1000000 

705 U1-2FT/T1 

710 U2-2FT/T2 

715 V3-2FVT3 


725 REM I RESET VIEWPORT 1 WIHDOW, PRINT OUT RESULIS I 

730 REM - 

735 UIEMPORT 0,130,0,100 
740 WIHDOH 0,130,0,100 
745 MOVE 73,98 

750 PRINT "ULTRASONIC HAMEFORM" 

755 MOME 70,43 

760 PRINT "CROSS-CORRELATION '{IF 

765 MOME 0,98 

770 PRINT "DATE! "|Dt 

775 PRINT "JJSPECINEN: *|SF 

?ee PRINT “THICKNESS- "|Tj" cn" 

785 PRINT "JTRAN9DUCER: "|TF 

790 PRIHT " JMETHOD - VELOCITY, DELAY," 
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793 

808 

803 

8)0 

813 

820 

825 

830 

835 

840 

843 

850 

855 

868 

865 

870 

875 

880 

885 

890 

895 

900 

903 

910 

915 

920 

925 

930 

935 

940 

943 

930 

953 

960 

965 


PRINT ’ cfi/uS uSJ* 

PRINT USING ^70, 20,10. 3D, 28, 10.30": -OUERL0P"," “,Ul f - %T1 
PRINT USING *7l?»20,ljD. 3D, 20* ID. 3D" : ’PHASE *,* ",V2," *,T2 

PRINT USING ■70,20,10.30,20, ID. 30*:*C0RR SVS," «,T3 

PRINT •JJ’IAll" POINTS /' WAVEFORM" 

PRINT USING “ 1,60 » 10. 40, 30* I * T I Ml NG ERROR ♦U_*, 03*1000000, * uS" 
PRINT USING "60, ID. 2D, 20": "OR ♦H_" » 1/C4C100, • 

REN 

REN I SHOW OPTIONS t MAKE DECISION ABOUT PROGRAM FLOW I 

REM — 

PRINT " JJJJJJJINPUT CHOICE SHOWN BELOW: J" 

PRINT ' fr TO REPEAT SAME SAMPLE " 

PRINT *'N' 10 START 0 HEW SAMPLE" 

PRINT " S' OR <ENTER> TO STOP * " 

PRINT ’ JCMOTCE 5 GC"1 
IHFOT 2l 

IF 2*=’R‘ THEN 215 
If Zfe-N" THEN 180 
END 

REN 

REN t SUBR0U1INE II — ON SRQ POLL CHECK I 

REM - 

POLL D, S| I 
RETURN 

REM - 

REM I SUBROUTINE *2 -- DRAW AXES J, W0UEFORH I 

REN 

CALL "MIN", 0, Ml, LI 
CALL 'MAX', 0, 112, L2 
Ml* INTO.GT < 83101 ) > 

W2-0! 

IF 0BS(M1KM2 THEN 9€5 
M2— Ml 
GO TO 978 
MI— M2 


970 H3*LGHM2-ML> 

975 W3-INT<M3> 

980 N4-INK<M2-HI>/10tW3)*l 
985 M3*-W4/'2 

990 WINDOW 0,W?,M3T10TW3, <W4+M3)*J0TU3 
995 AXIS W2/ie,W4*10tW3/-8,0,M3*J0tU3 
1000 2t*STR< 10TW3) 

1005 MOUE 1,CW4+M3)*10TW3 

101O PRINT •KHUHdbbBbHHHt!" 1 2*1 " "fU* 

1015 FOR 1-0 TO 8 STEP 2 

1020 MOVE l,W4tl0tW3^8tI*H3*l0rW3 

1825 PRINT USING " ""UytiUbbbb"',40.2D ,, :W4/8*l4M3 

1030 NEXT I 

1035 7*»STR< tOT.Mtt 1000000) 

1040 MOVE 01/1.35,M3*'10tW3*l. 1 
1045 PRINT •JJ"T7*;i" "JH*| 

1050 FOR 1-P TO I,0';STEP 2 
1055 MOVE I *H2> / 1.0, f13*,10tW3T 1 . 05 

1060 PRINT USING ' " "UbJ" " , 20. ID,S": l*W2/10*03/10tMl 

1065 HEXT I 

1070 CALL "DiSP’,0 

1075 RETURN 
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TABLE I. - COMPARISON OF DIGITAL METHODS FOR ULTRASONIC VELOCITY MEASUREMENT 


Method 

Metal sample 3 

Composite sample** 

Overlap 

Phase-slope 

Cross-correlation 

Time 

delay, T, 
us 

Random 0 
time delay 
error, 
±us 

Ultrasonic 
velocity, v, 
cm/ys 

Time 

delay, T, 
vs 

Random 0 
time delay 
error, 
±yS 

Ultrasonic 
velocity, v, 
cm/ys 

0.90873 

.90922 

.90873 

0.00012 

.00005 

do 

0.5500 

.5497 

.5500 

0.74850 

.72417 

.74975 




a Metal sample: 250-grade maraglng steel. Thickness: 2.499 mm. Estimated error In 

velocity: ±0.2 percent. Measured with a SO MHz center frequency transducer. 

^Composite sample: graphite/epoxy, 8-ply, {0 ± 45) symmetrical. Thickness: 1.125 mm. 

Estimated error In velocity: ±0.2 percent. Measured with a 5 MHz center frequency 

transducer. 

c T1me delay errors are based on standard deviation of time delay for 30 echo pairs at 
fixed location on material sample. 

d The cross-correlation method calculated the time delays to be Identical for all 30 pairs 
of echoes. 














/-BACK SURFACE 




(a) Cross-section showing principle echoes. (b) Time domain trace of principle echoes. 

Figure 1. - Specimen-transducer configuration for pulse-echo velocity measurements and principal echoes. Echo FS 
is returned by front surface of specimen and B1 and B2 are successive echoes returned by the back surface. Time. 
T. is the "round trip" time delay between echoes B1 and B2. 



(a) Schematic diagram for acquiring two successive echoes, B1 and B2. in the same time window. 
Figure 2. - Diagrams for two alternative waveform acquisition schemes. 



(bl Schematic diagram for acquiring two successive echoes, B1 and B2. in separate time windows. 
* Figure 2. - Concluded. 
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ECHO B1 ECHO B2 



h W H 

Figure 3. - Representative back echoes B1 and B2 separately windowed shown in digitized form. Time 
delay, w. between windows is a predetermined, known quantity. Time delay T, between echoes B1 
and B2 is calculated from T • W + (T2 - Tit by using one of the methods illustrated in figures 4. 5, 
and 6. 



Figure 4. - Result of digital overlap method for determin- 
ing delay between echoes B1 (solid) and B2 (dotted) using 
echoes shown in figure 3. 



VOLTAGE. mV 



FREQUENCY. MHz FREQUENCY. MHz 


(c) Amplitude spectra. Id) Phase spectra. 

Figure 5. - Typical back surface echoes. 81 (solid) and B2 (dotted), and 
their amplitude and phase spectra. In-window time delays based on the 
digital overlap method gives T1 ■ 39. 8 nsec and T2 ■ 91. 0 nsec. The 
phase-slope method gives T1 - 38. 9 nsec and T2 • 90. 1 nsec. Differences 
between the time delays can be attributed to the data sampling interval 
error of 0.4 nsec. 


ECHO B1 ECHO B2 



TIME, nsec 


(a) T1 • T2. 

Figure 6. - Result of generating the normalized cross-correlation function for three cases of echo 
positions. The displacement of the absolute maximum of the cross-correlation function relative to 
the zero reference gives the time delay C. For these ideal echoes C was equal to T2 - T1 as mea- 
sured by digital overlap within the data sampling Interval error of Q 4 nsec 
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T2 - T1 - 99. 2 nsec 
C - 99.6 nsec 
T- W + C 
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TIME, nsec 
<b)Tl<T2. 

Figures. -Continued. 


CROSS-CORRELATION FUNCTION 
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CROSS-CORRELATION FUNCTION 



(cl T1 > T2. 


Figures. -Concluded. 
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L50X10' 1 



(a) Time domain trace of echoes B1 and B2. (b) Normalized cross-correlation function 

T2-T1* 1748 ±1002 usee. for B1 versus BZ C- 1750 ±1002 usee. 



TIME, usee TIME, usee 


(c) Time domain trace of echoes B3 and B4 id) Normalized cross-correlation function 

T2-T1- 1766 ± 0.002 usee. for B3 versus BA C - 1 752 ± a 002 usee. 

Figure 7. - Examples of digital cross-correlation method for determining time delay between simultaneously 
windowed, successive echo pairs exhibiting increasingly low signal-to-noise ratios. Backsurface echoes 
for a graphite epoxy specimen using a 5 MHz transducer. Each pair of echoes is the average of 10 wave- 
forms. 
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